Conductance fluctuations in graphene devices with superconducting contacts in 

different charge density regimes 
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Conductions fluctuations (CF) are studied in single layer graphene devices with superconducting 
source and drain contacts made from aluminium. The CF are found to be enhanced by super- 
conductivity by a factor of 1.4 to 2. This (near) doubling of the CF indicates that the phase 
coherence length is 1$ > L/2. As compared to previous work, we find a relatively weak dependence 
of the CF on the gate voltage, and hence on the carrier density. We also demonstrate that whether 
the CF are larger or smaller at the charge neutrality point can be strongly dependent on the series 
resistance Rc, which needs to be subtracted. 



INTRODUCTION 

Graphene has proven to be interesting for the study 
of conductance fluctuations (CF) in disordered systems 
[1-5]. Some of its remarkable properties are the linear 
low energy band structure [6] and the K-K' degeneracy, 
which leads to unusual scattering [7]. In common sup- 
ported graphene devices disorder is inevitably introduced 
by the substrate and CF have been shown to be major 
corrections to the transport at low temperature. 

In numerical calculations [8, 9] stronger CF were found 
away from the charge neutrality point (CNP). Addition- 
ally, the CF were reported to be non-ergodic and they 
were suspected to result from variations in trajectories, 
rather than in phase shifts [8]. The results and interpre- 
tation of [8] were debated in [10], where CF were studied 
analytically and a similar magnitude of the CF was found 
but with fulfilled ergodicity. 

Superconducting electrodes can provide additional in- 
formation on the transport mechanisms in disordered sys- 
tems [11]. Phase sensitive Andrecv reflection can occur 
at the interface between the superconductor (S) and the 
normal conducting graphene (G). At energies below the 
superconducting gap A electrons from G can only enter 
S in the form of a Cooper pair. This charge transfer 
needs time reversal symmetry, which can be broken by 
the application of an external magnetic field. 

In the following, we look at CF in a two-terminal 
graphene device with superconducting contacts. We 
compare the CF in the superconducting and normal state 
of the electrodes. The latter is achieved by either increas- 
ing the temperature or by applying a magnetic field. The 
CF near and far away from the CNP are presented. 



SAMPLE PREPARATION AND 
EXPERIMENTAL SET-UP 

Mechanical Microcleavage [12] was used to deposit the 
graphene flakes onto a highly p-doped silicon wafer, cov- 
ered by 300 nm silicon dioxide. The silicon wafer serves 
back-gate where a voltage V ga te is applied. 



The graphene flakes on the substrate were located with 
optical microscopy and categorised by optical contrast. 
To ensure that the selected flakes are indeed single layer 
graphene, Raman spectroscopy [13] with 532 nm laser 
light was performed. 

Subsequently, a PMMA mask was used to pattern the 
electrical lines by electron beam lithography. A thin 
metal film (5 nm Ti / 50 nm Al) was evaporated on the 
sample in a UHV c-gun evaporation system. After lifting 
the PMMA mask off, the sample was cleaned in chemical 
solvents (acetone, ethanol and 2-propanol) and annealed 
in vacuum at 200 °C for several hours. 

The inset of figure 1 shows a micrograph of the device 
investigated. The single layer graphene flake is approxi- 
mately of length L = 1.3 /im and width W = 3.8 /im. 

Electrical measurements were performed in an Oxford 
Instruments 3 He cryostat at temperatures in the range 
of 230 mK - 1.5 K. A Stanford Research Systems SR830 
lock-in amplifier was used to apply a 20 /iV AC voltage 
which was superimposed onto a DC voltage V s d- The cur- 
rent through the sample was detected by an IV converter. 
No series resistors (like the electrical lines of the cryostat 
or the impedance of the IV converter) were subtracted 
from the data, unless mentioned. 



RESULTS AND DISCUSSION 

Figure 1 shows the measured two-terminal conduc- 
tance G through the device as a function of the back-gate 
voltage Vg a te, for two different temperatures, 230 mK and 
1.5 K. At low temperature, the Al contacts are in the 
superconducting state, while for the higher temperature 
they are normal. The charge neutrality point (CNP), 
indicated by the minimum in G, lies at V ga te = 15 V 
which suggests p-doping. Charged impurities introduce 
scattering and thus disorder [14] which is responsible for 
the broad shape of the dip in G. 

At large doping, i.e. for small negative gate voltages, 
the conductance G is reduced in the normal state by 
1 e 2 /h, while G remains the same near the CNP. This 
can be explained by the series resistance of the Al leads 
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FIG. 1. (online colour at: www.pss-b.com) The response 
of the conductance G through the single layer graphene de- 
vice on the applied back-gate voltage V ga te at 230 mK and 
1.5 K (lower curve). The charge neutrality region lies at 
Vgate = 15 V. The marked regions (I)- (III) are investigated 
later on in greater detail. The inset shows a micrograph of 
the device. 



adding to the device in the normal state. Figure 1 also 
shows that the CF are enhanced in the supercondcuting 
state. However, since the measurement in this figure 
was performed with an applied source-drain bias volt- 
age V s d = 140 /iV which is similar to A, the CF are not 
fully enhanced by Andreev reflection. 

We want to focus on three distinct regimes that can be 
identified in figure 1: (I) high doping with slow change in 
G, (II) intermediate doping with steeply changing G and 
(III) near the CNP. We expect diffusive transport as in 
disordered metals in regime (I), whereas in regime (III) 
near the CNP the graphene properties should dominate. 



near the CNP as a function of V s d and V ga te at 230 mK. 
In figure 2(a), the electrodes are in a superconducting 
state. Due to the superconducting gap A transport is 
expected to be dominated by Andreev reflection at low 
bias voltage V s d < 2A/e, where the factor 2 accounts 
for the two superconductors defining source and drain 
contacts. In this bias window the device displays a sup- 
pression G, indicative for weakly coupled contacts. The 
fluctuations in G as a function of V ga te are largest at 
zero bias and are reduced by an increased V s d- When 
applying a small magnetic field, as shown in figure 2(b), 
the superconductivity is suppressed. Consequently, the 
fluctuations of G at zero bias are visibly reduced. 

In order to quantify the change in the conductance fluc- 
tuations originating from enhancement due to Andreev 
reflections, we investigate the standard deviation of the 
conduction, SG as a function of the bias voltage V s d, as 
shown in [5] . 
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FIG. 3. (online colour at: www.pss-b.com) Conductance fluc- 
tuations SG as a function of the applied bias voltage V s d 
at 230 mK in region (III). Two regimes can be observed in 
the superconducting state (solid line): Outside of the super- 
conducting gap (V s d > 0.3 mV), the fluctuations are reduced 
in comparison to the values in the gap. In 100 mT perpen- 
dicular field superconductivity is suppressed and SG (dotted 
line). At V 3 d > 0.6 mV the normal and superconducting SG 
are identical. Below V s d = 0.3 mV, the superconducting con- 
tacts enhance the fluctuations by a factor of up to 1.5 when 
compared to SG in the normal state at 100 mT. 

Figure 3 shows a plot of 5G against V s d near the 
CNP in (III). When the electrodes arc in a super- 
conducting state, the CF are largest at zero V s d with 
SG = 0.22 c 2 /h. With increasing V s d the CF decrease 
and reach at V s d = 1 mV a value of SG = 0.14 e 2 /h. 
By comparing the slopes of the steeply changing CF at 
small V s d and the more slowly varying CF at larger V s d, 
we can extract a crossover and thus a superconducting 
gap of 2A re 300 [iV. In the normal conducting state at 
100 mT, the CF are about 0.14 e 2 /h at V sd = 1 mV, 
which coincides with the value of the superconducting 
electrodes. However, for V s d < 0.3 mV the CF in the 
normal state saturate at SG = 0.15 c 2 /h. Thus the CF 
in the superconducting state are enhanced by a factor of 
up to 1.5 compared to the normal state. 

When comparing the normal state fluctuations SG « 
Figure 2 displays grey scale plots of G in region (III) 0.15 c 2 /h to the theoretical Altshuler-Lee-Stone (ALS) 
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FIG. 2. (a) Grey scale of the conductance G as a function 
of the gate voltage V ga te and the applied source-drain bias 
V s d at 230 mK in region (III). The system enters the super- 
conducting gap for \V s d\ < 0.3 V. The conductance fluctu- 
ations in the gap region increase, while outside the fluctua- 
tions are decreasing with increased \V s d\ ■ A cut in the super- 
conducting gap (indicated by the dotted line) is shown by the 
black curve, (b) Measurement of G with a 100 mT perpen- 
dicular magnetic field applied. The conductance fluctuations 
are reduced, as can be seen from the cut (black solid curve) 
at V sd = V. 
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FIG. 4. (online colour at: www.pss-b.com) (a) The conduc- 
tance G as a function of the gate voltage V ga te at 230 mK 
with and without the subtraction of a constant series resis- 
tor Rc ■ When Rc is not subtracted (black curve) the CF are 
stronger in region (III) than in region (I). With Rc = 1.35 kfi 
subtracted the CF are stronger in region (I) and decrease 
when approaching the CNP in region (III), (b) Circuit dia- 
gram of the device. A constant resistor Rc is in series with 
the graphene resistor Rc- All fluctuations 8 R originate in the 
graphene part. 



value of SG = GMyW/L e 2 /h w 1.2 e 2 /h [8], wc find 
that the prediction is almost one order of magnitude 
larger than the measured value. One way to resolve this 
discrepancy is to assume that the phase coherence length 
lef, is much shorter than the length L of the device. As- 
suming that each segment of length 1$ is fluctuating in- 
dependently, yields an overall reduced CF [5]. Applying 
this approach to the current values yields a phase co- 
herence length of 300 nm. However, the fact that the 
measured CF is nearly doubled in the superconducting 
states suggests that almost all electrons are within the 
coherence length of one of the superconducting interface, 
hence 1$ > L/2. This in strong contradicting to the value 
of 1^ estimated from the normal state fluctuations. 

Another parameter which can reduce the apparent CF 
is a built-in series resistance Rc- This is schematically 
shown in 4(b). The graphene part has the resistance Rg 
which fluctuates with SR, whereas Rc is assumed to be 
constant. In practice Rc consists of the electrical lines 
in the measurement set-up, the input impedance of the 
IV converter and the resistance of the Al contact lines. 
In addition, it may also contain to some extent the con- 
tact resistance of the device [15]. All this is intrinsic to 
two-terminal measurements. We assume that Rc is in- 
dependent of the carrier concentration in the graphene 
flake. The subtraction of Rc corrects the measured fluc- 
tuations, 5G m , by 5G = (1 + Rc / Rg) 2 ■5G m . By adding 
the known resistance of the measurement setup (mea- 
surement lines and IV converter) we arrive at a lower 
boundary of Rc ~ 1-35 kf2. An upper boundary of 
Rc ~ 1.5 kO. can be extracted from the quantum Hall 
effect, which due to the magnetic field measures the de- 
vice in the normal state. Rc has been deduced from the 
shifts of the Hall plateaus. 

Figure 4(a) compares G against V ga te a t 230 mK with 



and without Rc = 1.35 kJ7 subtracted. As expected the 
correction results in the largest change at large doping 
when G is large. The same is true for the CF, as we will 
show in the following. 

In order to assess the sensitivity of the CF on the value 
of Rc, we have investigated the dependence of the CF 
measured at V s d = and 230 mK for different Rc values, 
ranging from to 2 kil in regime (I) (high doping), (II) 
(intermediate doping), and (III) (CNP). Figure 5 shows 
the result. In order to calculate SG a B-spline was sub- 
tracted to account for the non-constant background of 
the graphene conductance G(V gate )- The absolute values 
of the CF are lower than the values found in figure 3, 
because there only a linear background was removed. 
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FIG. 5. (online colour at: www.pss-b.com) Conductance fluc- 
tuations SG at 230 mK and V s d = in dependence on the 
gate Vgate and for different values of Rc- A B-spline was 
subtracted from the original data to account for the changing 
conductivity of the graphene. (a) Normal conducting elec- 
trodes at 100 mT. The effect of Rc is strongest in region (I). 
Without the subtraction of Rc, SG increases from (I) to (III). 
At the largest subtracted Rc = 2000 fi, SG decreases from 
(I) to (III), (b) SG in the superconducting state. The CF are 
enhanced compared to the normal state. The dependence of 
the CF on V ga te and Rc follows the same tendency as in the 
normal state. 

In figure 5(a) the electrical contacts are in a normal 
conducting state due to the application of a 100 mT mag- 
netic field. The as measured data (no Rc subtracted) 
shows the smallest CF which increase from 0.03 e 2 /h in 
(I) to 0.08 e 2 /h in (III). Around the series resistance that 
we estimated for our set-up (1.35 kf2 - 1.5 kfi), the fluc- 
tuations are nearly independent of V ga te, with values of 
~ 0.14 e 2 /h in regimes (I) and (III). If R c = 2 kfl is 
subtracted the CF in (I) reach 0.45 e 2 /h and are reduced 
to 0.2 e 2 /h in (III). 

Without the subtraction of Rc, the CF arc minimal 
at large doping and increase near the CNP. When Rc 
is subtracted the CF increase. However, they increase 
the most in region (I), less strong in region (II) and the 
weakest increase is in region (III). This can lead to a 
reversal of the dependence of the CFs on the V ga te- 

The conductance fluctuations 5G with super- 
conducting contacts are plotted in figure 5(b). The 
measured data without any Rc subtracted shows an 
increase from 0.07 c 2 /h in (I) to 0.16 e 2 /h near the 
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CNP in (III). Again, for values of Rc that lie within 
the estimates for our set-up, the CF vary only weakly 
with V ga t e from 0.32 e 2 /h in (I) to 0.25 e 2 /h in (III). At 
higher values of Rc the CF decrease significantly from 
(I) to (III). 

Evidently, the superconductivity enhances the CF in 
all three regimes. In (I) and (III) the enhancement lies 
close to a factor of 2. Yet in regime (II) it appears to be 
consistently higher and reaches a factor of 3. In part this 
large enhancement factor originates from the resistance of 
the aluminium electrodes, which is approximately 100 
in the normal state in our device. If we consider this 
difference, we find an enhancement factor of around 1.9 
for (I) and (III), but still above 2 for regime (II). 

In recent experiments contradicting observations were 
made. In [16] an increase of the CF near the CNP was 
seen, where the increase of the CF from high charge car- 
rier concentration to near the CNP was attributed to the 
formation of electron-hole puddles. In contrast, a de- 
crease of CF near the CNP was found in [17, 18]. The 
changes in the CF were reported to be of a factor of 3 or 
larger. Gate independent CF were observed in [19]. As 
in our device, [16] rely on two-terminal measurements, 
which include Rc, whereas [17, 18] use a four-terminal 
set-up, where Rc is irrelevant. As shown above, we can 
see all three gate dependences if we subtract Rc- Nev- 
ertheless, for the range of Rc that we estimate for our 
set-up, we find nearly no change in CF at high charge 
carrier density and near the CNP in our graphene de- 
vice. 

CONCLUSIONS 

Using graphene devices with superconducting contacts, 
we have shown that CFs are enhanced in the super- 
conducting relative to the normal state of the contacts. 
This enhancement amounts to a factor between 1.4 and 
2 and thus suggests that all electrons in the device of 
length L = 1.3 /mi are within the coherence length of 
one of the two contacts where Andreev reflection occurs 
which is the origin of the doubling [11]. The phase coher- 
ence length at 230 mK is therefore large and amounts to 

> L/2. We established a constant series resistor Rc 
to account for our two-terminal measurement set-up and 
device and investigated its impact on the CF. For values 
of Rc that we estimated for our set-up and device the CF 
remain smaller than e 2 /h. Furthermore, we investigated 
the dependence of the CF on the back-gate voltage V ga t e - 
A strong dependence on Rc is found. Without subtrac- 
tion of Rc the CF increase around the CNP compared 
to the CF at high doping. In the range of appropriate 
Rc the CF are nearly independent of V ga t e - 
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